*STATA Version 15.1

*Load figure1_replication_dataset.dta

recode v_56 v_120 (77=.)

ge income=v_120
recode income(1=2)
lab def inc_lb 2"<750" 3"750" 4"1000" 5"1250" 6"1500" 7"2000" 8"2500" 9"3000" 10"4000" 11"5000" 12"7500" 13">10000" 
lab val income inc_lb

ge taxsupport=v_56
recode taxsupport (1 2=1) (nonmissing=0)
logit taxsupport i.income i.c_0032
margins, over(income)
marginsplot, plotregion(margin(medlarge)) ysize(2.5) scheme(lean2) recast(scatter) ylabel(0(0.1)1, gmax labsize(medlarge) ///
             gstyle(dot)) xlabel(,labsize(medlarge)) yline(.5, ) title("") ytitle("% Support raise taxes for rich", size(large)) ///
			 xtitle("Net monthly household income, EUR", size(large))
